Size scaling of the addition spectra in silicon quantum dots 
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We investigate small artificial quantum dots obtained by geometrically controlled resistive con- 
finement in low mobility silicon-on-insulator nanowires. Addition spectra were recorded at low 
temperature for various dot areas fixed by lithography. We compare the standard deviation of the 
addition spectra with theory in the high electron concentration regime. We find that the standard 
deviation scales as the inverse area of the dot and its absolute value is comparable to the energy 
spacing of the one particle spectrum. 

PACS numbers: 73.23.-b, 73.23.Hk, 73.63.-b 



Measuring the current as a function of gate voltage 
provides decisive information on the addition spectrum of 
electrons in quantum dots. The addition energy is given 
by A 2 (A0 = E N+1 + E N -i - 2E N = ea5V g (N), where 
SV g (N) is the difference in gate voltage between the N th 

and the (N + l) th current resonance, a = ^ the ratio 
between the gate and total capacitance of the dot and 
En the ground-state energy of the A/"-electron state Q. 
Few experiments concent rate on the fl uctuations of the 
addition energy, a 2 = \J (A§j — (A2) 2 , which require a 
large number of resonances to be measured 0, EJ 0, H, . 
The distribution of A2 found in experiments is generally 
close to a Gaussian with a standard deviation cr 2 varying 
from values comparable to the mean one-particle level 
spacing Ai to much larger values. Until now the statis- 
tical variation of addition spectra were studied in dots 
made with high mobility 2D electron gases with a large 
kpl parameter, where J^f is the Fermi wavevector and i 
the elastic mean free path. The size of the dots varied by 
orders of magnitude between the different experiments. 
On the theoretical side, addition spectra have attracted 
considerable interest. Restricting ourselves to the case of 
a large number, A 7 ", of electrons, interactions beyond the 
constant charging energy model give a Gaussian distribu- 
tion of A2(Af) and increase the standard deviation of the 
addition spectra with respect to the constant charging 
energy approximation [H Q, H, E| 13 EI • These inter- 
actions are characterized by the ratio r s between the di- 
rect Coulomb and the kinetic Fermi energies. The result- 
ing standard deviation of the addition spectra includes a 
term due to the fluctuations of the single particle spectra 
and a term arising from the fluctuations of the charging 
energy Ec = e 2 /C^. The former contribution is esti- 
mated for chaotic or diffusive dots from the random ma- 
trix theory (RMT): a 2 ~ 0.52 Ai @. For small (r s < 1) 
or moderate values of r s 7], the latter is typically esti- 



* Electronic address: msanquer@cea.fr 



mated as 02 — Ai-^=, where gr oc ^- x {kpf) is the local 

dimensionless 2D conductance in the dot Jl [llJ| . As a re- 
sult (72 — Ai, even after introducing spin or exchange 
|l| effects in the calculations. Smaller values of gr result 
in a larger 02 12]. In the case of very large r s (negligi- 
ble kinetic energy), one expects a Maxwell distribution, 
where the standard deviation scales as the charging en- 
ergy (a 2 (xE c ) [HQ. 

We have analyzed the distribution of the addition ener- 
gies as a function of the area of the quantum dot, restrict- 
ing ourselves to the regime of large carrier densities in the 
dot. In contrast to previous experiments, we have stud- 
ied the standard deviation of the addition spectra in low 
mobility silicon quantum dots obtained by resistive con- 
finement from thin silicon-on-insulator (SOI) films [TEj . 
We find that the fluctuations of the addition energy scale 
as the inverse area of the dot. The fluctuations are com- 
parable to the mean level spacing for the one-particle 
spectrum, in agreement with most recent theories [Toj . 

The devices were made from boron doped (10 15 cm -3 ) 
SOI (100) wafers. In order to ensure a low resistivity for 
the electrical contact to the silicon, the silicon film lying 
under the contacts is thick (70nm) and heavily doped. 
During doping, the active areas are protected from amor- 
phization by a local thermal oxide layer. After its re- 
moval, we proceed to ion implantation (As ions above 
10 19 at. cm -3 ) of the thin active areas (12 to 22 nm thick). 
A hybrid deepUV/Ebeam lithography combined with re- 
sist trimming is used to pattern silicon nanowires with 
width W in the range 20 to 400 nm and length Lf. Af- 
ter the silicon wires have been etched, a 4 nm-thick gate 
oxide was thermally grown before chemical vapor deposi- 
tion of the in-situ doped poly-Si gate. Ebeam lithography 
was then used to pattern gates with length L down to 40 
nm. The back end sequence follows a standard CMOS 
process. Figure [3 shows a scanning electron micrograph 
(SEM) of the cross section of a typical wire. 
The confinement of the dot is obtained by the part of 
highly resistive SOI wire away from the gate. There is 
no deliberate oxidation of the silicon film to create tunnel 
barriers as in the PADOX process 16]. We used neither 
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TABLE I: Samples used in this study. W: wire width; L: gate 
length; Lf\ wire length; d: wire thickness; (SV g ): average 
peak spacing in gate voltage at T=4.2K; Ai: one-particle 
mean level spacing. The dimensions are in nm. 
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(SV g ) (mV) 


Ai (meV) 


SI 


70 


40 


200 


12 


3.8±0.3 


0.17 


S2 


50 


40 


200 


17 


5.9±0.4 


0.19 


S3 


70 


80 


200 


22 


2.8±0.1 


0.07 


S4 


50 
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22 


3.7±0.3 


0.09 
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0.13 
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FIG. 2: Peak spacing SV g versus gate voltage V g of a Si- wire 
transistor (sample S6) at T=4.2 K. The inset shows typi- 
cal conductance resonances as function of V g at T=l K and 
T=4.2 K. For the statistics only the high V g regime (many 
electrons) was considered, where 8V g is approximately con- 
stant (slope of linear fit (line): — 10 -4 ). 



FIG. 1: Left: SEM picture of a specific morphological charac- 
terization structure: cross-section of a 12 nm thick and 60 nm 
wide silicon nanowire. Right: Layout of a sample showing the 
2D electron gas (dash lines) of surface area A — WxL+2xdxL. 



fluctuations of the channel width to shape the quantum 

dot iiiiiiini, nor control gates to define the dot elec- 
trostatically. The latter procedure induces a deformation 
of the dot, which influences the addition spectra |flL l20j| . 
The area and shape of our dots is simply given by the 
overlap between the gate and the SOI film; the area can 
be as low as 3300 nm 2 and as high as 14400 nm 2 . This 
results in values for Ai, which are larger than in previous 
studies (see Table ID |1 S H . 

A 2D electron accumulation layer is formed near the 
surface of the nanowire (see Figure Q). The mobility was 
estimated from the sheet resistance of SOI films without 
gates. Under the gate it could be either enhanced, be- 
cause of the protection provided by the gate or decreased 
by electrostatic disorder at the gate oxide/poly-Si inter- 
face. From r\j c± 4000 Q (at room temperature as well as 
at T=4.2 K, for d=22nm) and the doping level, we esti- 
mate the mobility to /i ~ 150 cm _2 V _1 s _1 and (kpf) — 
3-5. r s — (aeohr Jirns ) ~ 1 is 0.6 for ns=10 13 cm -2 (1.7 
for ns=10 12 cm -2 ). 

Hundreds of samples with various geometries, doping 
and gate oxide thicknesses were integrated on 200 mm 
wafers. We chose 20 samples to study at low tempera- 
ture, which were cleaved and mounted in a He3/He4 di- 
lution refrigerator with a working temperature range of 
70 mK < T < 4.2 K. An ac drain-source voltage Vds was 
applied and the Ids current measured with a standard 



lock-in technique. Lossy microcoaxes were used for RF 
filtering. Due to the device variations (see Table 0}, we 
were able to adjust the access resistance to favor the ob- 
servation of single-electron transport. Under these con- 
ditions the resistances varied little between T=300 K and 
T=4.2 K with values of the order of 100 kft measured at 
large gate voltage. This value is typical for having resis- 
tive confinement as demonstrated in [2]| and analyzed in 
[22^ . The recorded features are remarkably stable, even 
after several cooling cycles. The statistics are obtained 
by single V g scans covering several hundreds of peaks. 

The inset of Figure [21 shows a detail of a typical drain- 
source conductance as a function of the gate voltage V g . 
The spacings SV g between the peaks for the whole gate 
voltage range are shown below the inset. Neglecting the 
first few dozen peaks, the mean value of the peak spac- 
ing in Figure [21 is (£Vg)=5.9±0.2 mV. An overview of 
(SVg) for different wire geometries is given in Table [J We 
have chosen samples with the same level of doping (10 19 
As/cm 3 ) and the same gate oxide thickness (t ox =4 nm). 
The slope of a linear fit in the range 0.5 < V g < 2 V is 
of the order of 10 -4 , indicating that the gate capacitance 
C g is little affected by the increasing number of electrons 
in the dot [2^. This is in contrast with the observation 
of Ref. , where a shift of the period with gate voltage 
was attributed to an electrostatic coupling between two 
dots. We compared the measured values of (SV g ) with 
the expected values SV ca \ = 77- for a planar capacitance 

C g = e r eoA/t ox , where A is the surface area of the 2D gas 
(see Figure 0} and e r is the relative dielectric constant of 
the oxide. Figure [3| shows this ratio between measured 
and calculated mean peak spacings for the samples listed 
in Table [I] The ratios are close to 1, which shows that 
the dot size is, indeed, close to the area A covered by 
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FIG. 3: Ratio of the measured (SV g ) and calculated (SV ca i) 
peak-to-peak distance for a flat capacitance C — eoe Y A/t x 
(e r =3.9, area A, oxide thickness £ x=4 nm). The inset shows 
measured oscillations after subtraction of the background in 
units of e 2 /h as a function of the gate voltage V g for three 
different areas (lowest curve: sample S6; middle curve: sample 
S4; highest curve: sample S3). Curves for S4 and S3 were 
shifted for clarity. 
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FIG. 4: Fluctuations of the peak spacing distributions (72, 
normalized to the mean one-particle energy level spacing Ai, 
versus area A (filled dots: 4.2 K, open dots: 1 K). The inset 
shows the typical Vds — V g diamonds of the first conductance 
peak in a A=3300 nm 2 sample. The one-particle levels appear 
as additional lines in the diamonds. Projecting along the 
V g axis and taking a = 0.4 gives the measured excitation 
spectrum, with a mean value Ai = 0.6 ±0.2 meV, in good 
agreement with the calculated values of Table U] (see text). 



the gate. We estimate the mean level spacing Ai = 2 ^ A 
(m*=0.19m e is the effective mass of the electron), taking 
into account the spin and band degeneracy for electrons 
on the (100) silicon surface |24j |. We compared this value 
with direct measurements of the excitation spectra in a 
sample of comparable size (A = 3300 nm 2 ) but different 
doping from samples SI to S7. For low carrier densities 
we find additional levels in the typical diamond shaped 
Vds — Vg dependence shown in the inset of Figure 21 
The measurement of the level distances is in good quan- 
titative agreement with the estimation of Ai considering 
that, due to the so-called Lifschitz tail, Ai is expected to 
be larger at lower densities, by roughly a factor of two or 
three 13. 

For the analysis of the distribution of A2 we assumed 
a constant value for a, as in previous works, although it 
decreases slightly with V g . This variation is enhanced at 
low V g (in the first peaks region omitted from the statis- 
tics) because of the reduction of the source and drain 
capacitances. As usual in MOS-Field Effect Transistors 
near threshold, the capacitances between the channel and 
source/drain and body, respectively, increase with the 
gate voltage, while the gate capacitance remains constant 
|2q . a is obtained from the slopes of the Coulomb di- 
amonds at intermediate gate voltages and by fitting the 
temperature dependence of the resonances at high gate 
voltage, a c± 0.4 is found for all samples from the analysis 
of the Coulomb diamonds. We also note that a can fluc- 
tuate from resonance to resonance around its mean value, 
an effect neglected here which deserves a more system- 
atic study. Figure [5| shows the distribution of the peak 
spacings in gate voltage SV g for three different areas. As 



shown in the inset of Figure 03 the distributions shift to 
higher mean values of SV g as the sample size decreases. 
At the same time a broadening of the distribution can be 
observed. The measurements were usually done at T=4.2 
K, i.e. ksT > Ai. A comparison with measurements 
at lower temperature, T=l K and T=300 mK (effective 
minimal electronic temperature), has shown that the ad- 
dition spectra is sligthly modified: 02 increases at lower 
temperature, but this increase lies within the error bars of 
the T=4.2 K statistics (see inset of Figure |5J). As in pre- 
vious reports we observe a Gaussian distribution for the 
addition energies with a variance 02 , which is comparable 
to Ai. The Gaussian shape contradicts the Constant In- 
teraction approximation within the RMT model, where 
the ground state energy is the sum of the charging energy 
and the energies of the occupied single particle states [l| . 
To compare our results with theory [l(|, we plotted 02 
normalized to Ai as a function of the dots' area (see 
Figure 2]). Note that the simulation in Ref. 10] uses pa- 
rameters comparable to ours: rs — 1.5, gr — 4 and 
N ~ 200. Due to the smaller size of our dots, Ai is order 
of magnitudes larger than in @, 0- Alternatively, the 
inset of Figure [51 shows the standard deviation of the 
frequently used normalization (SV g — (SV g )) / (SV g ) , ob- 
tained in Gaussian fits. The standard deviation increases 
when scaling down the area of the dots. reaches about 
10 % for the smallest sample (A=3300 nm 2 ). This value 
is close to the ratio Ai ~0.19 meV over the value of the 
addition energy, A 2 = ea{6V g )=2A meV: Ai/A 2 =0.08. 
The normalization of 02 with Ai (Figure leads to a 
constant -g- (here —1), which proves that 02 scales with 
Ai in our samples. At low gate voltages it is obvious from 
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FIG. 5: Relative peak spacing distribution P(SV g ) for three 
different dots' areas at 4.2 K. The solid lines are the gaussian 
fits used to extract the standard deviation. The shift towards 
larger values of SV g for decreasing areas is clearly visible. The 
distribution becomes larger when the sample area decreases. 
This enhancement of the fluctuations in smaller samples is 
shown in the inset, where the standard deviation a of the 
normalized fluctuation (8V g — (SV g ))/(SV g ) is plotted against 
the surface area of the sample. 



the localized regime. As the added electrons and their 
screening clouds are both localized on different sites, dis- 
tant from almost zero to the diameter of the dot, the rel- 
ative fluctuations of the addition spectra can reach 100% 
as explained in Ref.|l3q. 

We presented a systematic study of the size depen- 
dence of the addition spectra in low mobility and ultra 
small silicon dots. Our dots are very similar to ultimate 
silicon MOSFETs, different from those used in previous 
systematic experiments. The small size implicates a one- 
particle mean level spacing of the order of a few Kelvin, 
which is much larger than in previous experiments. Due 
to the chosen fabrication technique, the dot shape is not 
modified by the gate voltage, which only modifies the 
concentration of electrons in the dot. We have found an 
excellent agreement with theories in the same range of 
parameters, both for the expected values and the size 
dependency of the fluctuations. The presented analysis 
is restricted to the high concentration range of carriers. 
Experiments investigating the low density case, where 
the dots undergo the metal-insulator transition, are in 
progress. 
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